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Abstract 

We construct new gravity backgrounds holographic dual to neutral plasma with 
U(l) global symmetry in the presence of constant electric field, considering its full 
back-reactions to the metric. As the electric field and the induced current cause a 
net energy in-flow to the system, the plasma is continually heated up and the cor- 
responding gravity solution has an expanding horizon. After proposing a consistent 
late-time expansion scheme, we present analytic solutions in the scheme up to next- 
leading order, and our solutions are new time-dependent solutions of 5D asymptotic 
AdS Einstein-Maxwell(-Chern-Simons) theory. To extract dual CFT stress tensor 
and U(l) current from the solutions, we perform a rigorous holographic renormaliza- 
tion of Einstein-Maxwell-Chern-Simons theory including full back-reactions, which 
can in itself be an interesting addition to literatures. As by-products, we obtain 
interesting modifications of energy-momentum/ current Ward identities due to the 
U(l) symmetry and its triangle anomaly. 
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1 Introduction 



AdS/CFT correspondence provides an important opportunity for both gravity and strongly 
coupled CFT [T]. Many features in gravity that have been found before the correspondence 
find their nice physical re-interpretations as field theory phenomena in the dual CFT. 
Study of a strongly coupled CFT, especially its time-dependent dynamics, in terms of 
conventional field theory techniques, is hampered by many intrinsic difficulties. AdS / CFT 
correspondence is one consistent framework that can provide reliable information on some 
aspects of these strongly coupled gauge theories. 

Finite temperature system, or plasma, seems particularly interesting to study in the 
framework [21 [3], because of its possible relevance to quark-gluon plasma one observes 
in heavy- ion collisions which is believed to be strongly coupled The expectation is 
that some finite temperature properties might depend weakly on details of microscopic 
theories, and be universal as long as interactions are strong [5], [6]. In the gravity side, a 
finite temperature plasma is conjecturally described by a black-hole, a space-time with 
horizon. Dynamics of the horizon seems to describe long wave-length hydrodynamics of 
the dual CFT plasma [21 [TJ El |9], while theory-specific non-hydrodynamic phenomena 
would involve the whole non-linear bulk space-time dynamics. There has been much 
recent interest in these subjects; see Refs. [IDl ttH [12] for reviews. 

One way of studying plasma is to perturb the system and observe its linear responses. 
This gives quite a lot of useful information on the properties of the plasma, mostly related 
one way or the other to hydrodynamic properties. AdS/CFT correspondence is equipped 
with formalisms which nicely fits to this particular need . They are boundary pertur- 
bation [H] and holographic renormalization [15]. As the perturbations one imposes are 
small and treated linearly, one typically forgoes considering back-reactions to the metric. 

A less studied aspect is going beyond linearized approximations and taking back- 
reactions to the gravity equations of motion. One typically needs this extension in order 
to study far-from-equilibrium dynamics of the plasma, such as early-time thermalization 
or black- hole formations [TH [T71 [HI [191 [201 [211 [22]. It is also conceptually more appealing 
to have a full complete solution of the system, rather than working in an approximation. 
Another reason for including back-reactions is to study effects to the energy-momentum 
tensor from other symmetries/operators of the theory. For example, external gauge po- 
tential coupled to a global symmetry of a CFT can give rise to violation of conformal 
Ward identities. To rigorously compute such effects, one needs to consider back-reactions 
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to the metric. 

In this work, we initiate a study for one particular case that indeed needs full back- 
reactions to the metric. Our situation is a CFT plasma with U(l) global symmetry 
in the presence of constant external electric field. Our model includes N=4 SYM with 
U(l) R-symmetry as a particular case. Because CFT possesses excitations of arbitrary 
small energy scale, the applied electric field inevitably causes an induced current along 
its direction. Previous linearized study of this situation gives us a conductivity [23] (see 
Ref. [21] for one non-linear study), but our intention is to go beyond this simple linearized 
description. Observe that the electric field and the induced current do work on the system 
and there is a net energy in-flow to the system, 

±^E.J . (1.1) 

SO that the plasma will be continually heated-up. To correctly describe what is happening 
in a more long-time basis, one necessarily has to take into account gravity back-reactions 
from the U(l) sector. Especially, a natural expectation is to have an expanding hori- 
zon corresponding to heating the plasma, which might be an interesting time-dependent 
gravity solution in itself. 

The action we study is 

{IQnG,) C = R+ 12- -FmnF^'"" - ^—e'^^^^^AMFr^pFQ^ , (1.2) 

where k = for N=4 SYM. Although it is not difficult to write down simplified equa- 
tions of motion of the situation exploring relevant symmetries of the system, such as 
translation-invariance in the present work, it seems in general to be a hard task to find 
exact analytic solutions. Instead, we follow a less ambitious path by invoking a consis- 
tent late-time expansion scheme, so that one can solve easier sets of differential equations 
order by order. The idea is similar to the Janik-Peschanki's late-time expansion of boost- 
invariant plasma [25], although black-hole horizon in our case is expanding rather than 
receding as in their case. We stress that identifying correct scaling variable as well as the 
right expansion parameter that we do here for our case is not a trivial result. Within the 
proposed scheme, we succeed in finding smooth analytic solutions up to next-leading order 
in the expansion. Although the results will be more and more complicated as one goes 
beyond higher orders, there is no further conceptual block in our scheme. Our solutions 
seem new to our knowledge. 

Because our solutions are taking full back-reactions to the metric, one also needs 
a full holographic renormalization [TB] of the Einstein-Maxwell(-Chern-Simons) system 
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including back-reactions. To our knowledge, there has not been complete analysis of 
this, and we intend to fill this gap in this work. Particularly interesting by-products 
of our analysis are a set of modifications of conformal Ward identity as well as energy- 
momentum/current conservation laws. From these results, we finally obtain the CFT 
stress tensor and the current of our new solutions up to next-leading order. 

It is worth of mentioning that a similar problem with a constant magnetic field, instead 
of an electric field, was studied by D'Hoker-Kraus some time ago [26], with a quite different 
motivation. A major difference in their case is the static nature of the solution, which 
simplifies analysis in a fundamental way. In order to be more complete, we also tried 
to search for a static solution with constant electric field without current to observe a 
singularity at a finite radius, which makes the solution unacceptable. We think this is 
physically understandable because the current will be inevitably induced and the energy 
in-flow must happen so that the geometry is necessarily time-dependent l]. 

We can broadly divide our paper into three parts. In the first part (section 2 and 
section 3), we propose a consistent late time expansion scheme and subsequently solve 
the equations of motion in the particular scheme till the first leading order. The next 
part (section 4) of the paper is devoted to using a rigorous holographic renormalization 
approach to derive the modified trace anomaly and modified Ward identity in the presence 
of an external electric field and taking into account the full backreaction. This can be an 
interesting result on its own independent of the other results of the paper. In the third 
part (section 5), we use the solution derived in section 3 to extract the energy momentum 
tensor and current density and observe that the obey the modified Ward identity derived 
in section 4. We don't think our analysis is complete in any respect, but our intention 
is to report our current findings which might be interesting to a further study on the 
subject. We finish our paper by giving a few open, future problems to extend the current 
analysis. 



2 Constructing gravity solution : late time expansion 

The situation we are aiming at is a time-dependent finite temperature plasma with a U{1) 
global symmetry of which a constant external electric field is applied to the system. For 
simplicity, we will be considering a neutral plasma, and we leave the case of finite charge 
density as a future direction. We are interested in obtaining an explicit 5-dimensional 

-'^We thank Andy O'Bannon for raising this question to us. 
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time- dependent gravity solution that is holographic dual to the situation. Note that 
time- dependence is an unavoidable feature in the presence of the electric field, because a 
current J is expected to arise in response to the electric field and there will be a net 
energy flow into the system proportional to J ■ E, by which the plasma will be heated up. 
Consequently, one expects the dual gravity solution to have an increasing event horizon 
as time goes on. 

As one expects the plasma at hand to be heated by the applied electric field, it is 
natural to have a thermodynamic concept of temperature, at least in sufficiently late time 
regime. The problem of finding an exact gravity solution with a constant electric field is 
an interesting subject itself, although it seems difficult to our eyes to handle analytically. 
Instead, we intend to study a late time asymptotic solution invoking a systematic late 
time expansion that will be discussed shortly, and to solve the equations of motion order 
by order. We have checked that our expansion scheme works consistently at a few lowest 
orders, and it seems very likely that the scheme is consistent at all orders. The philosophy 
is similar to the case of late time expansion of boost invariant plasmas [25], while the 
difference is that boost invariant evolution is a cooling process and the event horizon is 
receding, contrary to our case of heating up. 

To set the notation, let a homogeneous electric field point to x^-direction with a 
constant magnitude E. This will fix the unique boundary condition of the bulk U{1) gauge 
field at the UV boundary of 5-dimensional asymptotic AdS space we are constructing. 
One advantage of having a neutral system is that while the electric field will induce 
a net current J = J^x^ along the direction of electric field, there won't be any net 
momentum flow generated along the field. Intuitively this can be understood by the fact 
that a neutral plasma has an equal amount of thermally excited, positively and negatively 
charged carriers that will move in opposite ways under an external electric field, so that 
there will be no net momentum flow generated by the current carriers. More rigorously, a 
neutral plasma has a charge conjugation symmetry under which the electric field and the 
current flips its sign while the momentum flow is left unchanged. Invariance under this 
charge conjugation requires the momentum flow to vanish. The absence of momentum flow 
will greatly simplify our gravity solution ansatz: one can safely put the {t, x*) [i = 1, 2, 3) 
components of the metric to zero. Using diffeomorphisms one can further choose the 
gauge where 

9rr = gri = , 1 = 1,2,3 , Qtr = 1 , (2.3) 

where r is the 5'th holographic coordinate. Geometrically it is an affine coordinate of 
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ingoing null geodesies in Eddington-Finkelstein coordinate. After renaming (x^, , x^) = 
{x,y,z), the resulting gravity ansatz that is consistent with a left-over 5*0(2) symmetry 
of [x^jx"^) rotation is written as 

ds^ = -r^a{r,t)dt'^ + 2drdt + r%{r,t)dz'^ + r^c{r,t)(^dx^ + dy^^ , (2.4) 

with three unknown, yet to be determined functions (a, b, c) which depend only on r and 
the time t. We adopt the Eddington-Finkelstein coordinate as it has been proven to 
be useful to implement the correct causal boundary condition at the horizon by simple 
regularity of the metric components [7] . It was also shown that late time expansion of boost 
invariant plasmas works consistently at all orders in this coordinate [27], and we conjecture 
the same in our case too. To have an asymptotic AdS space with flat 4-dimensional CFT 
metric at the boundary, one needs to impose the UV boundary condition at r — )■ oo 

(a,6,c)^ (1,1,1) , r^oo , (2.5) 

for all t. 

Regarding an ansatz for the bulk gauge field, it is clear that one needs to turn on 
component only, in the gauge At = 0. The non-normalizable mode of encodes the 
external electric field by the UV boundary condition 

A.^Et , r oo , (2.6) 

while its normalizable mode, which should be determined by solving equations of motion, 
contains the induced current Jz along the field direction. Note that one will have to 
turn on At as well in the case of charged plasma, which would complicate the analysis. 
Therefore, the gauge field of our interest can be put simply as 

A = A{r,t)dz , (2.7) 

with one more function A{r, t) whose UV boundary condition is 

A{r,t) ^ Et , r^oo . (2.8) 

With the above Ansatzs for the metric and the U{1) gauge field, the equations of motion 
become a set of partial differential equations of four unknown functions (a, b, c, A) with 
the specified boundary conditions. With a lack of further intuition to these differential 
equations in finding an exact solution, we will instead invoke a systematic late time 
expansion scheme that allows us to solve them order by order analytically. 
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To identify the correct scaling variable and expansion parameter, we first observe that 
the applied electric field becomes more perturbative to the plasma as it gets hotter as 
time goes on, because due to conformal symmetry what matters is an effective dimen- 
sionless strength of the electric field, where T is the temperature. As T increases 
one therefore expects that adiabatic linear response theory will describe the system more 
accurately. Sub-leading effects coming from non-linear responses will be suppressed at 
late time regime by powers of an expansion parameter (^^^ where the square is due to 
the charge conjugation symmetry. In linear response theory of an electric field applied to 
an adiabatic plasma of temperature T, the induced current is proportional to the electric 
field and is given by 

J, = a{T)E , (2.9) 

where ct(T) is the electrical conductivity at temperature T of a neutral plasma. Conformal 
invariance dictates that <j{T) oc T. Then the energy in-flow due to the applied electric 
field will be 

^r^T^^^J-Er^E^-T , (2.10) 

dt dt ^ ^ 

where we have used the fact that a conformal plasma energy density e has ~ behavior. 

The easy solution of the above differential equation for T is 

r~^i-P , (2.11) 

for the leading late time behavior of the plasma at large t. In the gravity side, temperature 
T roughly corresponds to the location of black-hole horizon, that is, rn ~ T, and in the 
leading dual gravity solution one may expect to have a horizon moving towards the UV 
region as rn ^ t^ in the late time regime. This naturally motivates a scaling variable 

u=^ , (2.12) 

in terms of which the leading solution will look "static", and there will be sub-leading 
corrections to this leading solution by suitable negative powers of t. A natural candidate 
for this sub-leading expansion parameter is the previously mentioned non-linearity 
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yj^2j ' (2-13) 

and we conjecture that this late time expansion scheme works consistently order by order. 
Notice also that the temperature at late time grows more and more slowly 

dT 2 , , 

— ^ra^O , (2.14) 

Cit' 
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so that our adiabatic analysis in the above seems self-consistent. 

More explicitly, our first naive proposal, which in fact needs a refinement that will be 
discussed shortly, would be 

{a,b,c){r,t) = J2i(^n{u),bn{u),Cn{u))t~^ , u=— , (2.15) 

n>0 t3 

for the functions entering our metric ansatz with the boundary conditions as m — > oo, 

(ao(M),6o(M),Co(M)) ^ (1,1, 1) , (a„(ix),6„(«),c„(«)) ^ (0,0,0) , n>l .(2.16) 

One needs a slightly more care for A{r, t) in the gauge field. The leading expectation for 
the current Jz is 

J, a{T)E ■ E r^t^ , (2.17) 

and this current sits in the normalizable (~ ^) mode of Az, so that the leading late time 
behavior of A{r, t) near UV boundary should be 

n 1 

A{r,t)r^Et + — + --- = Et+^ + --- . (2.18) 
This suggests the following late time expansion for the function A{r, t), 

A{r,t) = Et + -Y.jMt-'r , (2.19) 

n>0 

with the boundary condition, 

jn{u) ^0 , n>0 , u^oo . (2.20) 

The validity of the expansion scheme f l2.15p and fl2.19p can be tested by performing a 
few lowest order terms in the equations of motion to see whether one finds a consistent 
solution order by order. Although one finds a good solution at leading zero'th order 
which is presented explicitly in the next section, there seems to appear an inconsistency 
at subsequent sub-leading orders, which necessitates an improvement of our first guess 
above. One however finds that the equations of motion indeed have the right expansion 
parameter fl2.13p . so that our arguments for (tt,t~t) seem more or less correct up to a 
possible subtlety yet to be identified. 

One observes a hint for the resolution in a near boundary analysis of general solutions 
of our equations of motion, which one typically performs in the procedure of holographic 
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renormalization [15]. We will discuss more details about holographic renormalization in 
section m Expanding (a, b, c, A){r, t) near r — > oo with the boundary conditions (12. 5 p and 
(12. 8p . one finds that the equations of motion dictate 



a(i)(f) (a«(t))' 1 6(4) (^) ^2iQg^ 



6(r,t) = 1 + 



, 2 



r 4 r"^ Y 2 24 y 12 

A(r,t) = Et + :^+f/^)(t)-ii?a«(t))l + 0-i^ + ---, (2.21) 



with four undetermined functions {a^^\ a^^\ b^^\j^'^^) of time with one constraint equation 

d,a^'\t) = j^'\t) , (2.22) 

which will be shown to lead to energy conservation law. The function a^^\t) is in fact 
spurious as it arises from a simple coordinate reparameterization 

r^r + ^a^^\t) , (2.23) 

and subsequently one safely ignores it from now. {a^^\b^^\ j^'^'^) encode the expectation 
values of energy-momentum and U{1) current in the dual CFT whose precise formulae will 
be given in section [5] later, but what is interesting for us at the moment is the appearance 
of several terms in the metric functions (a, b, c) due to a backreaction of the external 
electric field E. In the case of pure Einstein gravity, this term is given by certain 
curvature tensors of boundary CFT metric which is related to the conformal anomaly, 
and in our case of fiat CFT metric it would have vanished. In the presence of external 
gauge potential that couples to U{1) current, it gets additional contributions from field 
strengths of the external gauge field as our results indicate. We stress that the presence 
of terms is necessary for a consistent solution of equations of motion. Moreover, as 
the equations of motion are non-linear, one generically expects higher powers like 

/ log r ' 



j,4 



m>l , (2.24) 



to appear in further sub-leading expansions of actual solutions. 

What we observe in the above paragraphs, (which we will discuss in more details in 
our section on holographic renormalization), is that in the holographic renormalization 
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procedure, presence of log terms in the near boundary solutions is inevitable in the pres- 
ence of non trivial electric field and non trivial background CFT metric. The late time 
expansion scheme is supposed to approximate the full bulk solution (which includes the 
near boundary behaviour as well) at late times and hence it should commensurate well 
with this observation of near boundary behaviour of the fields. This calls for refining the 
late time ansatz by addition of suitable log terms. In the next few paragraphs we will 
identify such terms. 

As discussed in the previous paragraphs, in terms of our scaling and expansion vari- 
ables {u = rt~^,t~t), a consistent expansion scheme of (a,6, c, A) should be able to 

4 

reproduce these logarithmic structures near r — )■ oo. In terms of (u,t~3), one writes 
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logr logii 4 1 , , , 

and one observes that the first piece can easily be captured by next leading terms of the 
previous naive expansion scheme, while the second piece with logt can never be taken 
care of because there is no way to introduce logt in the naive expansion scheme. The 
easiest way to remedy the problem is to invoke a refined expansion scheme including 
powers of (\ogt ■ in addition to usual powers of We conjecture and confirm at a 
few lowest orders that logt indeed comes in a combination (\ogt ■ t~^^, so that multiple 
logs are suppressed by at least same number of powers of t~3. This naturally leads to our 
proposal of late time expansion scheme 

(a,&,c)(r,t) = J2 («(n,m)(«),6(n,™)(^),C(„,„)(ii)) rt("+™) (logt)" ,u = ^. (2.26) 



(n,m)> (0,0) 



One arrives at a similar conclusion for A{r,t), too. From the equations of motion, one 
finds that there appears a logarithmic term in the near boundary analysis starting at 
whose coefficient is roughly of order of dtj^'^\t). As we know that j^'^\t) ~ (see 
eq. fl2.18p ). we have a term of 
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2 log r If log M 4 1 , , , ^ 
-[^t-^ + —\ogt-t-^] , 2.27 

which guides us to the following expansion 

A{r,t) = Et + \ J2 Jin,m)iu)t-I^^+"^^ {hgtr . (2.28) 

(n,m)> (0,0) 

At each order > 0, all functions of total order n + m = N will enter the equations 
of motion in general. Note that while the number of unknown functions increases as A^, 
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there are also number of equations to solve because the equations of (logt)™' factors 
with m < N are all independent. Eqs. (12.261) and fl2.28p are our proposal for a consistent 
late time expansion of the gravity background corresponding to heating a plasma by an 
external electric field. 

We finish this section by commenting one caveat. The equations of motion have an 
obvious time-translation symmetry of t — )■ {t + to), so that from any solution one gets 
another by simply replacing t with {t + to). In particular our scaling variables {u,t~^) 
maps to 

(u,r^^) ^ (u- (l + tot-^) ' ,t-l ■ (l + tot-') , (2.29) 

under this transformation. Starting from our expansion (I2.26P and (12. 28 p . this transfor- 
mation would generate a new series of sub-leading terms with additional integer powers 
of t~^, so that a more general series that also covers this freedom seems to be 

^ rl("+'»)rp(iogt)'" . (2.30) 

(n,m,p)>(0,0,0) 

However, it is natural to conjecture that these additional terms with are rigidly deter- 
mined with one free integration constant t^, and reorganize themselves to our expansion 
of (I2.26P and (I2.28P with t being (f + to) instead, so that one can always return to our 
proposal (I2.26P and (I2.28P by suitable translation of time coordinate t. We leave a proof 
of this claim to the future. 

3 Leading and next leading order gravity solution 

In this section, we provide a consistent analytic solution in our late time expansion scheme 
introduced in the previous section up to the order n + m < 1. We expand the Einstein 
equation = and the Maxwell equation Ai^ = in powers of (t~^,logt), and 
let us denote the coefficient of t~"(logt)'' by a superscript {a, [3). Then £^^^'^\ S*^u'^\ 
gu{i,o)^ gu{o,o)^ gz{o,o)^ gx{o,o) ^ g^j^^ A^i^'°^ provide complete differential equations for the 

leading order functions of (n, m) = (0, 0), and one also checks that these equations are not 
affected by sub- leading functions of (n, m) > (0,0). As these leading order equations are 
non-linear, there is no systematic way to solve them, and one generally needs an educated 
guess. Once the leading order functions are found however, the differential equations for 
sub-leading functions in further expansions of equations of motion are linear differential 
equations, and one can solve them order by order systematically. 
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Because one expects the leading solution to be a black-hole solution whose event 
horizon would look static in the scaling coordinate u, one can try the following ansatz 
guided by the static black-hole solution, 



i(0,0) 



u 



1-^ 



Vo)(«) = C(o,o)(m) 



(3.31) 



with some constant w to be determined later, and one indeed finds that this solves almost 
all the above differential equations, namely, S^^f'^\ S*^u''^\ £^^z'^\ are 

satisfied by this ansatz. The remaining two equations £^^^'^^ and A^i^'°^ become then a 
first /second order differential equation for the remaining j(o,o)(''^) respectively. Explicitly, 
f gives us 



■^(0,0)1 



u 



while Aiz^' provides a second order differential equation 



u 



w 



•^(0,0)' 



u 



+ (Su' + 



w 



3(0,0) 



'u) + Eu^ 







(3.32) 



(3.33) 



It is simple to integrate f l3.32p . but notice that the integrand on the right-hand side has 
a simple pole at the event horizon u = w which would cause a logarithmic singularity 
of j{o,o), unless this pole is removed by a suitable choice of constant w such that the 
numerator also vanishes a.t u = w. This uniquely fixes w to be 



w 



and the integration of fl3.32p with the UV boundary condition j(o,o)(oo) 



J{o,o)\ 



u = — 




(!)' + 



u 



(f) 



+ w 



(3.34) 
results in 

(3.35) 



It is then nicely checked that this solution also solves the remaining equation fl3.33p . 
which is a non-trivial test for consistency of the solution. In summary, horizon regularity 
and UV boundary condition uniquely determine the leading order solution without any 
further integration constant. This feature will in fact be true in subsequent sub- leading 
higher-order solutions too. 

The differential equations for next leading order functions of (n, m) = (1,0) or (0, 1) 
are provided by the same set of Einstein equations and Maxwell equation as above with 
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a — )■ a + I and /3 = or 1. As they are linear differential equations, solving them has no 

conceptual problem while it is algebraically quite cumbersome to present, so we only sketch 

/ t(- 0) t(- V 

the procedure and simply present the results. One first observes that iSu' ■,£ u' ' 

are second order differential equations for the combinations (&(i,o) + 2c(i o), ^'{o,i) + 2c(o,i) 
respectively, whose integration gives us 

/ 



^(1,0) + 2c(i,o) 



+ 



^(0,1) 



+ 2c, 



(0,1) 




— u\+-\og 



(f)'+ 



u 



(f)'+ 



\ 



(3.36) 



where we already fixed one integration constant by UV boundary condition while the 
remaining integration constants (Ci,Ci) will be discussed shortly. Next one finds that 
the equations (£^^1^'°'' + 2£^^i^'°'*, + 2£^x^'^^) provide first order differential equations 



of (a(i,o), 0(0,1)) respectively after using the above (13. 36 p . It is not difficult to integrate 
them to have 



«(1,0) 



«{0,1) 



+ 




log (^t + 2 



3U' 



(3.37) 



with two more integration constants (C2, C2). Then, with the above results of (I3.36P and 
(I3.37p . one can check that the equations (£^*j^'°\ f^*^/'^"*) and (£^"i^'°'*, f^"!^'"*^"*) are automati- 
cally satisfied. Next using (I3.36P and (I3.37p . the equations (£^'^1^'°'', £^^1^'^"*) become second 
order differential equations for (c(i,o), C(o,i)) respectively, whose first integration produces 



C(l,0)l^J 



+ 



3^^ 2>v? 




u 



8 
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with integration constants (C^jC^,) to be determined shortly. Note that the right-hand 

2 

sides of (I3.38P have simple poles at the horizon u = w = (^^^ unless (C3, C3) are chosen 
to remove the singularity. This regularity condition determines (C3, C3) uniquely as 



C3= (^yj (l + 21og(2i?3jj , C3 = . (3.39) 

If one wishes, one can further integrate f l3.38p with UV boundary condition to find 
(c(i,o), C(o,i)) explicitly. For C(o,i), the result is simply 

C(o,i) = ^ , (3.40) 

while the result for C(i^o) looks quite complicated and we skip writing it down explicitly 
by presenting only its near boundary asymptotic as 

C E'^ f 1 1 \ 

C(i,o)~^ + ^(31ogn-21og^+-log2--j+--- , w^oo , (3.41) 

for a later convenience. Combined with this, f l3.36p and f l3.37p give us complete next 
leading-order solutions of (a, 6, c) of (n, m) = (1,0) and (0,1) with integration constants 
(C'1,2, C'1,2) yet to be determined. 

The remaining final equations to solve, namely (£"|^'°'', £^"|^'^'') and (A^i^'^'', A^i^'^'') 
provide first/second order differential equations for (i(i,o) 5^(0,1)) respectively. The first 
order equations (£^"j^'°\ f^'^j^'"^^) take the following structure 




•^(1,0)' 



^(1,0) («) , U'-(l) )4i)(") = ^(o>i)W ' (3-42) 



where (F(i 0), F(o,i)) are some complicated combinations of {a,b,c) and j(o,o) so that they 

depend on the integration constants (Ci,2, (^1,2)- To avoid logarithmic singularity at the 

2 

horizon u = w = ("f)'') one has to require -F(i,o),(o,i)('"^) = 0- By explicit computa- 
tions, one finds that -F(i,o),(o,i)('"') do not involve (Ci, Ci) so that this regularity uniquely 
determines (6*2,(72) only. Explicitly one has 

eI , „ 4, „ 12 „^\\ „ . . Co 



^(1,0) («^) = rrr 1 - log 2 - 4 (log 2) ^ - - log + — (C2 - 3^2) , F(o,i) («;, 

323 V 3 E"^ ^ ' J 2sEs 



. 2 „ 1 ) 
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so that this fixes (C2, C2) uniquely as 

After fixing ((72,(72), it is not difficult to integrate fl3.42p to obtain j(o,i) as 

~ E u 
i{o,i){u) = -Ci-- j^^g j-^ , (3.44) 



while the expression for j{i,o) is too complicated to present here. For our later purpose, 
we present explicitly only the near UV behavior of 

Jii,o){u) ~ (1 - 2Ci) + O . (3.45) 

It is not difficult to check also that these results for j(i.o) and j(o,i) from f l3.42p satisfy the 
second order differential equations {Aiz^' ,J^z^' ), which is a consistent check for the 
framework. 

Finally, the remaining integration constants {Ci, Ci) can easily be identified as spuri- 
ous modes coming from a simple coordinate re-parameterization 



, Ci , Cilogt 
r ^r + — + — — — , (3.46) 
bt DC 



or equivalently 



Ci , Cilogt 

u u-\ J H 4 — , (3.47) 

which is a left-over gauge freedom in our coordinate system. See a^^\t) in eq. fl2.2ip . 
Therefore one can safely set them zero, and one is left with the unique next leading 
solution presented above with all the integration constants completely determined by 
regularity at the horizon. 



4 Holographic renormalization of Einstein- Max well- 
Chern-Simons theory 

The purpose of this section is to perform a rigorous holographic renormalization of 
Einstein-Maxwell-Chern-Simons theory, including full back reaction to the metric. Our 
motivation is two-fold. First, to our knowledge there has been no complete analysis on 
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a holographic renormahzation of Einstein-Maxwell theory taking into account back re- 
actions to the metric, and it seems useful to clear the situation including Chern-Simons 
term too. Our second purpose is to extract interesting physical observables from our new 
gravity background in the previous section, such as energy-momentum and U{1) current, 
using the rigorous results of holographic renormahzation. The necessity of a careful holo- 
graphic renormahzation is due to the presence of external electric field which breaks the 
underlying conformal symmetry, and one naturally expects a violation of conformal Ward 
identity T^^ = 0. This hints to a possible non-trivial effect to the energy- momentum 
tensor from the electric field, and at least conservatively, one needs to follow rigorous 
steps of holographic renormahzation including back reactions to the metric to correctly 
identify it. 

As we will formulate regularization and counter-terms in a covariant way, we can 
choose to work in the standard (and simplest) Fefferman-Graham coordinate, although 
our gravity background is written in an Eddington-Finkelstein like coordinate. After 
identifying covariant expressions of counter-terms, one can safely apply the results to 
our Eddington-Finkelstein coordinate (or any other reasonable coordinate) to extract 
renormalized finite values of physical observables. 

4.1 Near boundary solution to the equations of motion 

A la Ref.[T5], one starts with near UV boundary behaviors of the metric and U{1) gauge 
field which are dictated by the equations of motion, 

ds"^ = GMNdx'^^ dx^ = -\ — Qnyix, p)dx'^dx^ , 

p 

A, = 40) + 42)^ + i?(2)plogp+... , (4.48) 

with 

9,A^, P) = gf} + 9f}p + g^^p" + h^^p" log p + ■ ■ ■ , (4.49) 

where and A^^'^ are Dirichlet boundary conditions that are given by hand as external 
CFT data such as background 4-dimensional metric and an external potential coupling to 
f/(l) current. As in the case of pure Einstein gravity and Maxwell theory without back 
reaction, the terms g^'^\ h^'^^ and -B'-^^ are expected to be completely fixed by equations 
of motion in terms of these external CFT data, g^^^ and A^^\ The contributions from 
A^^^ especially to g^'^^ or h^^^ would be a new aspect coming from back reactions, whose 
presence can already be seen in our previous near boundary behaviors f l2.2ip . The next 
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terms of g^'^^ and A^'^^ are dynamical, and not completely determined by near boundary 
equations of motion, but are constrained by them. Finite, renormalized CFT observables 
are given in terms of these coefficients, and the constraint equations from equations of 
motion typically give us Ward-identities for these observables. As mentioned before, we 
are interested in these Ward-identities too. 

The action for the Einstein-Maxwell-CS theory is 



S 



WttG 



jj'xdp^ {r[G] + 12 - - ^e*™^AMF;vpFQ«' 



+T7rV / d'x^2K (4.50) 
IbTrGs JdM 

Where K is the trace of the second fundamental form (external curvature) and 7 is the 
induced metric on the boundary. The Gibbons-Hawking boundary term in the second 
line is necessary for the consistency of Dirichlet problem for manifolds with boundary. In 
the coordinate system fl4.48p the equations of motion read as, 

- \Tr{g-'g") + -Tr{g-'g' g-'g') + ^ [-ApTr{g-' A'^ A!) - -TT{g''Fg-'F)^ = 0, 

Wg'^,-V,Trig-'g') - pg'^'^A'^F,^ = 0, 



4p 



mc[g] + p [-2g" + 2g'g-'g' - Tr{g-'g')g') + 2g' + Tr{g-'g')g 
+ ^g{8p'A'^g-'A' - pTrig-'Fg-'F)) + ^pFg''F -2p'A'^A' = 0, 

^ d,{^g^''A',) - ^^e^'^-'F^^F^s = 0, 

= 0. 
(4.51) 



_ -9 ^ ' ^V-9 

■d, (V^g'-'K) - -^d^ {V^ig-'Fg-'r) - -^e^-^A'^F. 



Differentiation with respect to p is denoted by prime, is the covariant derivative 
constructed from the metric g^u{x, p) and Ric[g] is the Ricci tensor of the metric g^j.,^ 
We will now invoke the near boundary expansions as in ( I4.48P in the above equations 
f l4.5ip . The solution we obtain are 



^The convention for Riemann tensor that we foUow is R^Jpq — dpTji^Q + r^pF^g — P ^ Q which is 
opposite to that of [15] . Therefore there is a change in sign at various places in the equations of motion 
and in the analysis thereafter. 

•^From now on we will suppress explicit mention of in the terms. For example we will write 
F^^ = (Fg-^F)^^ and Tr{F^) = Tr{g-^Fg-^)F and so on. 
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^(2) 



(4.52) 



96 96 32 / 32 ^ (o)^^^!' ^ 3 ' 

rK^?(4)) = ^rr(^f2)) - ^rr(Ffo)), 

V^^it^ = V^(if/(V-^(rr^(2))^il) + ^((Tr^7(2))^-rr(^^^^ 
-V^ (^(Tr(Ffo)))e^) + Ig^'^'^ + FW. 

The last line in the above equation can be solved as 

= ^^?(2W - |(rr^(2))^?j + ^ ((7^r<?(2))^ - Tr{gf,^)) g^^J - ^{Tr{F^oM;j + V , 

(4.53) 

where t^,^ is arbitrary, but according to (14.52 p . satisfies the following constraint 

Trt = l(^rr(^(22))-(rr^?(2))^ + ^(TrFj))) . (4.54) 

4.2 Divergences and counter terms in terms of induced metric 

The on-shell gravitational action diverges near the boundary and hence we have to regulate 
the theory and remove potential divergence terms by adding local counter terms in the 
boundary. We follow the steps outlined in [15] . We first regulate the theory by restricting 
the bulk integral to the region p > e. The regulated action is given by 



S, 



1 



reg 



1 



+ 



IGttG^ JdM 



f d^xdpV^ (r[G] + 12 - \fmnF 
Jm \ 4 



'^MNPQRa p p 



d X\ 
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1 



1 



p3 12 Je p-" 



+ 



d X 



1 

'p2 



^ 1 d^x^^) (6e-2 + a(4) log e) + 0(e°) 



(4.55) 



Where a 



(4) 



■ {Trg(2)Y '~ ^ i^^(-^{o)) • -'■'^ order to get a finite action we have 

to subtract the divergences by adding local counterterms in the boundary. For that we 
have to invert the relations between the induced metric 7^,^ = \g^v and g^^} perturbatively 
in e. The results are 115] 



-^(0) 
Tr{g(2)) 



(l-^eTr(<,(,)) + 0(e(^))) , 

Using the above results the counter term action takes the following form 

Rb] 



Tri 



(2); 



IdttGr J 



6 H h (34 [7] loge 



(4.56) 



(4.57) 



lOvrGs 

Where a^b] same expression as given below f l4.55p with 5f(o) replaced by 7. While 

writing down the counter terms, we did not write down terms involving Rb]"^ ot (log e)-R[7]^ 
since these would yield finite contribution which would just be a different choice of scheme 
and will not affect physical quantitities like trace anomaly. 



4.3 Holographic stress energy tensor, current, and Ward iden- 
tities 

Having found the counter terms required to cancel the divergences of the action in the 
previous section, we are in a position to derive the expressions for the expectation value 
of the stress energy tensor and current in the dual theory which are given as 



-2 6Sr 



1 ^SfQYi 



lim(^-T,.[7] 



V-9io) SA 



(4.58) 



where 7)tiy[7] and J'^b] expressions for the stress energy tensor and current of the 

theory described at p = e without taking the limit e — )■ 0. From the above expressions it is 
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obvious that the expansions of ^^^[7] up to C(e) and J^[t] up to C(e^) will be important. 
Both of these quantities have contributions from Sreg as well as S^*. The expressions 
coming from Sreg are 
1 



"reg 



1 

^2 



8nG. 



deg^iu{,x, e) + g^y{x, e)Tr[g ^{x, e)deg{x, e)] - -gf,u{x, e] 



(4. 



0(e^e'log 



The contribution coming from the counterterms fl4.57l) are 



rpct 
f 

'J/ 



1 



fll/ 

Ct 



SttGs 
1 



V.F^^M loge 



where 



1 6 (/ d'^Xy/—ja4^['j]) 



rpa 



-7 



V.F'^'^[(?(o)]loge, 



2/^?J[7] • 



(4.60) 



(4.61) 



In order to apply the above results to fl4.58p . we need to expand them in powers of e, for 
which the following result will be usefulp 



1. 



1 



(0) 



6 



(4.62) 

After using this expression and after a slight bit of algebra, we observe that the ^ pole 
and logarithmic divergence in (Tf^u) and the logarithmic divergence in ( J^) cancel and we 
get perfectly finite expressions for both of them which follows 



1 



2V + 3/^it^ + ^Q(TrFj))^(°)-F(°) 



(4.63) 



The expression for (T^j^) obtained above is almost the same as obtained in [15] with the 
exception of the last two terms and explicit appearance of the boundary gauge field in the 
expression of t^^ and /i^fj. The trace anomaly that we obtain from the above expression 
is 



Tr{T) 



1 



1(4) 



167rG5 



(Tr^( 



2). 



(4.64) 



^The results are the sanie that appear in |15j with shght change in signs of some terms because of the 
opposite convention used by us. 
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Note that the terms proportional to /i^'^J and ^{TrF^QS^)gj^J — Fj^J"^ in (T^u) are scheme 
dependent and can be removed by addition of local finite counterterms such as 

j d^x^ [aa(4)[7] + /?Tr(Fj))] , (4.65) 

by appropriate choice of a and [3. Indeed as expected they do not contribute to the trace 
anomaly. The trace anomaly now contains a term depending on the boundary gauge field 
and even in the fiat space limit g^^^ — ?■ ?], we will still have a nonzero trace anomaly as 
expected. Removing the scheme dependent terms from the expectation value of the stress 
energy tensor and using the solutions to the equations of motion f l4.52p and f l4.54p . we get 
the following Ward identities 

V,(J^) = ^^=6^-^FWFi? . (4.66) 



Note the interesting contributions from the Chern-Simons term. 



5 Energy-momentum tensor and current of our solu- 
tion 

After developing the necessary tools in the previous section, required to derive the energy 
momentum tensor and U{1) current in the boundary from a bulk solution, we can apply it 
to our new gravity solutions up to next-leading order to find energy-momentum tensor and 
the current of our interest. Some results in the previous section are covariant, especially 
the counter-terms fl4.57p . so that they can be used in our Eddington-Finkelstein like 
coordinate in sections 2 and 3 as well. In our situation, the boundary CFT metric g'*^''-* 
is fiat, so various pieces such as ^f*-^-* are simply absent, which makes things a lot simpler 
than it looks. For completeness, we will state each steps of holographic renormalization 
in our coordinate system explicitly, while main necessary computations can be simply 
borrowed from the previous section. 

One first regularizes the action by considering only the bulk region of r < (note 
that p in the previous section corresponds to p = in our coordinate). One can easily 
check that p = e boundary in Fefferman-Graham coordinate is asymptotically identical 
to r = boundary in our Eddington-Finkelstein coordinate, so that the covariant 
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counter-terms in the previous section can be used without any modification. The exphcit 
expression for the counter-terms is given by fl4.57p . 

S"' = - / d'x (e + ^i?[7] + a4[7] log e) , (5.67) 
from which the counter-term contribution to the energy-momentum tensor is 

^^'^ = " 8^ i^^'" ' \ (^^^^^^ " \^^^^^^') + 2^?-^ log ' (5-68) 
with h^^"^ now becomes, using the fiat CFT metric g'^'^'> = rj, 

m = ^ {f^^^''''-') + ^F(°)F(°)% , (5.69) 

where 7 is the induced metric on the boundary, and A^^'^ is the external gauge potential 
coupled to the U(l) current. In f l5.69p . indices are raised/lowered by g^'^^ = rj. The 
regularized bulk contribution to the energy-momentum tensor is given by the Brown- 
York expression even in the presence of Maxwell field, [f| 

= ^ (^M. - ^7m.) , (5.70) 

which is already covariant so that we can also use this in our coordinate system. It 
is then straightforward to compute finite, renormalized energy-momentum tensor in our 
coordinate system using the above expressions, 

7;T = liml(T;:^(e)+T;t(6)) . (5.71) 

After plugging in the general, spatially-homogeneous near-boundary solutions of f l2.2ip 
in the presence of external electric field that we reproduce below 



a(r,t) = l + ^+Q(a«(t))'-a,a«(t)) 



1 ^ a(^)(t) _ EHogr ^ 

^2 ,y^4 ^ 



, a«(t) (a^^Kt)Y I 6(4) (t) EHogr 
b(r,t) = 1 + — + - — 7; + M ^ + 



, 2 



c(rt) ^ l + ^+(^l+f-V)(t)-^U + ^i^ + ... 
^'^ ^+ r + 4 r^^[ 2 24jr4+12 + ' 

Air,t) = Et + ^+U'\t)-^-Ea^'\t))\ + 0.^-^ + ..., (5.72) 
r V 2 / r"^ 



^This is because using the bulk equations of motion, the variation of bulk action with respect to metric 
variation always reduces to the same surface boundary term which is the Brown- York tensor. Maxwell-CS 
term does not produce any boundary surface term with respect to metric variation. 
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one obtains indeed the finite, well-defined CFT energy-momentum tensor as 



SttGs V 2 ' ' 6 

rpxx^j.yy ^ _J_(_}_ai^)u\i,Wft\_]_EA , (5.73) 
SttGs V 2 ^ ^ ^ ^ 24 / ' ^ ' 

with dynamically unknown functions a^^^ and b^^^ that can be read off from our explicit 
solutions in section 3. As claimed before, a^^\t) never appears in the results because it 
is a spurious mode of coordinate re-parametrization. Note that the above result nicely 
confirms the conformal anomaly (14 .64^ we obtain in the previous section. 



Tl^ = - R2 ^ I p(0) p{0)fiu (r 



which is a consistency check. Due to this constraint, one can't completely ignore contact 
terms of E"^ in the energy-momentum tensor above. However, as long as it satisfies the 
trace anomaly (15 .74^ . one is allowed to move them from place to place, which are different 
scheme choices. A similar procedure can be followed to obtain the renormalized value of 
U(l) current; note that the same counter term (14.571) also takes care of divergences from 
the Maxwell(-CS) action. The result is simply 

J. = ^3^'\t) , (5.75) 

where another dynamical function appearing in (I5.72p . From the constraint (I2.22p 

dictated by the equations of motion 

dta^'\t) = \e j^'\t) , (5.76) 

one also checks the Ward identity (I4.66P to hold, 

d.T^"" = F(°)^V, , (5.77) 

where the Chern-Simons piece in (I4.66P is not relevant in our specific situation. 

What remains to obtain final results of {T^^, J^) is to simply retrieve the necessary 
dynamical functions {a^^\h^^\ j^'^'^) from our explicit solutions in section 3, and it is not 
difficult to get 

a(^)(t) = -(^|)%t + ^(^logt + 21ogS + 6(log2)'-ilog2-2)+0(rt) , 
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/'\t) = -(|)*iUlt-+o(r^) . (5.78) 

These are the main results in this paper. In particular we see that the T** component of the 
energy momentum tensor increases with time |^ which is consistent with our expectation of 
energy being pumped into the system by an external electric field and the rate of increase 
of energy density being pumped into the system is also consistent with the Ward identity 
arising out of a rigorous holographic renormalization presented in section H] taking into 
account the full backreaction of the gravity solution. We also see that the energy density 
T** has the correct (T being the temperature) or t^^^ behaviour in the leading order as 
expected from a conformal plasma but what is non trivial are the subleading terms which 
are the new results of our analysis. 



6 Discussion and future direction 



Time dependent gravity background is a valuable play ground to check several interesting 
dynamical phenomena in gravity. It becomes even more interesting for the asymptotic 
AdS gravity as it holographically describes a dual conformal field theory at strong coupling 
whose study is typically intractable by conventional field theory techniques. Lattice stud- 
ies for time-dependent dynamical aspects also possess intrinsic problems. In this paper, 
we provide one more example of "quasi-analytic" time-dependent AdS gravity solution by 
applying external electric field to the system of U(l) symmetry, which pumps in energy 
flow into the system so that the plasma can be heated up continuously. It is "quasi- 
analytic" because the solution is obtained by invoking a late-time expansion scheme we 
propose, from which one can construct solutions order by order. We check consistency of 
the scheme and provide analytic solutions up to next-leading order, from which we ob- 
tain late-time behaviors of energy-momentum tensor and U(l) current. For this purpose 
we also perform rigorous holographic renormalization of 5D asymptotic AdS Einstein- 
Maxwell(-Chern-Simons) theory considering full back-reactions to the metric. Note that 
back-reactions should be included to find precise modifications of energy-momentum Ward 
identities due to U(l) symmetry. We find our results for the current density and energy 
momentum tensor to be consistent with the Ward identity derived in section HI We also 

^Note that the solution presented is rehable only in the late time regimes as we have worked out the 
solution in a late time expansion scheme (j2.26l2.28p 
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obtain the correct leading behaviour ~ T'^ for the energy density and also obtain non 
trivial subleading behaviours which are the new results in this paper. We would also like 
to stress that the results of section 4, on the modification of conformal Ward identities 
and trace anomaly in the presence of a background gauge field, are new and interesting 
on its own irrespective of the other results of the paper. Specifically we performed a 
rigorous holographic renormalization procedure on the Einstein-Maxwell-CS system by 
taking into account full backreactions from the metric and the results (14.641 14.66P gives 
us the complete modified conformal Ward identity and trace anomaly in the presence of a 
non trivial background gauge field. The appearance of the coefficient of the Chern Simons 
term k in the above results seems to encompass the effects from triangle anomaly in the 
field theory side and it would be interesting to study this effect by taking into account 
background magnetic field. We will discuss on this open issue more towards the end. 

We close this section by a few questions regarding future directions. 

Apparent horizon and entropy density : It is not the area of event horizon, but that 
of apparent horizon which gives one the entropy density in the dual gravity plasma [28] . 
This distinction is important in a time- dependent background. It will be interesting to 
compute entropy density as a function of time by identifying apparent horizon of our new 
geometry. 

Hydrodynamic interpretation with external electric field : Because the ratio ^ becomes 
arbitrary small in the late time regime, one expects that our dual plasma may be described 
by some form of hydrodynamics with electric field. Especially interesting question is 
whether our next-leading solution gives some kinds of transport coefficients. However, 
hydrodynamics is normally based on spatial-gradient expansion, while our plasma simply 
doesn't have any spatial-gradients because we are assuming translation invariance. This 
might imply that no hydrodynamic modes are excited at all, and non-hydrodynamic 
modes might be governing the dynamics of our plasma. The situation seems similar to 



Boost-invariant symmetry and early-time thermalization study : Instead of assuming 
translation symmetry that we study in this paper, one can also impose boost-invariance 
on the system with external electric field. It simply corresponds to a different prepara- 
tion/initial condition of the system. Boost-invariance frame naturally causes a cooling 
tendency due to metric expansion while electric field will try to heat up the plasma, so 

^We thank Andy O'Bannon for a discussion on this. 
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there can be interesting interplay between the two effects. Another advantage of hav- 
ing boost-invariance rather than translation-invariance is that a boost-invariant plasma 
is guaranteed to have hydrodynamic modes excited, so the question of hydrodynamic 
description can be more transparent [21]. Independently, it is also interesting to study 
early-time thermalization process starting at zero temperature and applying electric field 
after t = 0. These will be addressed in a near future [29] . 

Charged plasma with electric field : Although we restrict our focus on a neutral plasma 
in this work, it will be an interesting extension to consider a charged plasma with external 
electric field. The additional complication will be that the electric field causes the net 
momentum flow along its direction, and one can not remove (^tj-components in the 5D 
metric ansatz in any frame. Even 4D boost-transformations cannot completely remove 
them from the bulk 5D metric because charge conjugation symmetry no longer prohibits 
them. 

Both electric and magnetic fields, and effects from triangle anomaly : The case with 
both electric and magnetic field pointing to a same direction seems particularly interesting 
in regard to effects from triangle anomaly [301 Ell Ell ESI EH ESI ESI EZ] ; there will be a 
constant anomalous creation of charge density due to electro-magnetic fields, in addition 
to current/momentum flows. As the 5-dimensional Chern-Simons term is holographic 
dual to the triangle anomaly, it will play an essential role in this case. On a practical 
side, what one needs is a more general ansatz in the gravity side including {At,gu) as 
well, while Fu = B will be a constant in full 5D due to Bianchi identity. It seems to be 
a tractable problem to our eyes. 
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